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We investigate the electronic structure of cobalt atoms on a copper surface and in a copper host by 
combining density functional calculations with a numerically exact continuous-time quantum Monte 
Carlo treatment of the five-orbital impurity problem. In both cases we find low energy resonances in 
the density of states of all five Co d-orbitals. The corresponding self-energies indicate the formation 
of a Fermi liquid state at low temperatures. Our calculations yield the characteristic energy scale - 
the Kondo temperature - for both systems in good agreement with experiments. We quantify the 
charge fluctuations in both geometries and suggest that Co in Cu must be described by an Anderson 
impurity model rather than by a model assuming frozen impurity valency at low energies. We show 
that fluctuations of the orbital degrees of freedom are crucial for explaining the Kondo temperatures 
obtained in our calculations and and measured in experiments. 

PACS numbers: 71.27.+a, 75.20.Hr, 73.20.At 
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I. INTRODUCTION 

The Kondo effect, arising when localized spins inter- 
act with a metallic environment, is a classic many body 
problem^ At present, idealized models dealing with, for 
instance, a single spin degree of freedom screened by a 
sea of conduction electrons are well understood. Spin 
S = 1/2 Kondo models or single orbital Anderson impu- 
rity models have been widely considered to describe mag- 
netic impurities with open d-shells in metallic environ- 
ments and proved helpful in qualitative discussions. 2 ~ 9 
As realized, however, by Nozieres and Blandin in 1980, 10 
such idealized models may ignore important aspects of 
the nature of transition metal impurities as they disre- 
gard orbital degrees of freedom. This makes compar- 
isons between theory and experiment often very difficult. 
More realistic models accounting for the orbital struc- 
ture, Hund's rule coupling, non-spherical crystal fields 
and an energy- and orbital-dependent hybridization of 
the impurity electrons with the surrounding metal are 
theoretically very demanding due to the multiple degrees 
of freedom and multiple energy scales involved. 

In recent years different attempts have been made 
to address this problem. For the classic example of 
Fe in Au, which has been experimentally studied since 
the 1930s, a model describing the low energy physics 
has been derived 11 by comparing numerical renormal- 
ization group (NRG) calculations to electron transport 
experiments. The Kondo temperature Tk, below which 
the impurity spin becomes screened and a Fermi liquid 
develops, served as a fitting parameter in this study. 



A scaling analysis of multiple Hund's coupled spins in 
a metallic environment showed that Hund's rule cou- 
pling can strongly quench the formation of Kondo singlet 
states^ 2 - For highly symmetric systems like Co adatoms 
on graphenei^ or Co-benzene sandwich molecules in con- 
tact to metallic leadsp^ the orbital degree of freedom has 
been suggested to control Kondo physics down to the low- 
est energy scale. However, a general strategy to assess 
which degrees of freedom are involved in the formation 
of low energy Fermi liquids around magnetic impurities 
in metals is still lacking. 

Co atoms coupled to Cu hosts present another exper- 
imentally extensively studied system which has been in- 
terpreted in terms of Kondo physics £d£r£J£. Theoret- 
ical descriptions of this system have often been based 
on single orbital Anderson impurity models^— or Kondo 
models- and the role of orbital fluctuations in these sys- 
tems has remained rather unclear. Recently developed 
continuous time Quantum Monte Carlo (CTQMC) 16 ap- 
proaches allow to describe the full orbital structure of 
magnetic impurities in metallic hosts, while accounting 
for all electron correlations in a numerically exact way. 
So far, however, such CTQMC studies have been limited 
to rather high temperatures^ well above typical Kondo 
scales on the order of 10K to 500K. The realistic descrip- 
tion of transition metal Kondo systems thus remains a 
long standing open problem in solid state physics. 

Here, we employ the recently developed Krylov 
CTQMC method 18 in combination with density func- 
tional based first-principles calculations to achieve an ab- 
initio description of two archetypical Kondo systems: Co 
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Kondo temperature from our Quantum Monte Carlo data 
PVC|) . In Section [V] we discuss the prevalence of charge 
fluctuations and the implications of their presence (|V A[l . 
We elaborate on the role of spin and orbital fluctuations 
and their influence on the low energy behavior (|VB[) . 
Section |VI1 is a summary and conclusion emphasizing the 
implications of our findings and outlining prospects for 
future investigations. 



FIG. 1: (Color online) Sketch of Co impurities in a Cu host. 
We consider the two cases: (i) the impurity buried in the bulk 
(Co in Cu), and (ii) on top of the Cu layer (Co on Cu). 



adatoms on a Cu (111) surface, as well as Co impuri- 
ties in bulk Cu (see Fig. [T]). We consider the energy 
dependent hybridization of the impurities with the sur- 
rounding host material as well as the full local Coulomb 
interaction and find low energy resonances developing in 
the spectral function as the temperature is lowered. Such 
resonances are found in all impurity 3c?-orbitals and our 
calculations indicate that spin- and orbital-fluctuations 
are crucial for the formation of low energy Fermi liquids 
involving all impurity 3d-orbitals. We also demonstrate 
the intermediate- valence character of the Co impurity in 
bulk, which implies that the physics cannot be correctly 
described by a low-energy Kondo model which neglects 
charge fluctuations. 

The paper is organized as follows. Section [TT] defines 
the model we use to study the Kondo physics in Co 
coupled to Cu hosts. Section [TTT] specifies the Density 
Functional Theory calculations of the hybridization func- 
tion piIAp and describes the Krylov CTQMC method 
used to solve the five-orbital impurity problem (IIIIBI) . 
In Section IIVI we present the calculated quasi-particle 
spectra pV A[) . exemplify the low energy Fermi-liquid be- 
havior based on the self-energies pVB[) and estimate the 



II. MODEL 

A realistic description of the Co atoms on Cu (111) 
and in bulk Cu including all five Co 3d orbitals can be 
formulated in terms of a multi-orbital Anderson impurity 
model: 
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It describes an impurity characterized by quantum 
numbers a (orbital and spin), with corresponding an- 
nihilation operators d ai on-site energies e Q , and local 
Coulomb interactions U ai ,...,cn- This impurity is embed- 
ded in a sea of conduction electrons described by annihi- 
lation operators c k and dispersions e k , where k includes 
crystal momentum, band index and spin. The coupling 
between the impurity and the conduction electrons is pro- 
vided by the hybridization Vk a - 

In the AIM only the impurity site is subject to a quar- 
tic interaction term, whereas the bath of conduction elec- 
trons is assumed to be non-interacting. The bath degrees 
of freedom can thus be integrated out and the local elec- 
tronic properties of the impurity can be described by the 
effective action 
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with 

GaXxctz ( iuJ n) = {uo + n) - A QlQ2 (iu n ) (4) 

and the hybridization function 



performed first principles density functional theory cal- 
culations. 
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To specify the parameters of impurity models describ- 
ing Co on Cu (111) as well as Co in bulk Cu we have 



III. COMPUTATIONAL METHODS 

A. Density functional calculations 

Density functional theory (DFT) calculations were per- 
formed to obtain relaxed geometries and the hybridiza- 
tion functions for single Co atoms in Cu and on a Cu 
(111) surface. The DFT calculations have been carried 
out using a generalized gradient approximation (GGA)i£ 
as implemented in the Vienna Ab-initio Simulation Pack- 
age (VASP)2£ with projector augmented waves basis sets 
(PAW) i 21 ' 22 For the simulation of a cobalt impurity in 
bulk Cu we employed a C0CU63 supercell structure. Co 
on Cu (111) was modeled using a 3 x 4 supercell of a Cu 
(111) surface with a thickness of 5 atomic layers and a 
Co adatom on the surface, see Fig. [TJ All structures were 
relaxed until the forces acting on each atom were below 
0.02 eVA" 1 . For Co in bulk Co the entire supercell and 
for Co on Cu (111) the adatom and the three topmost Cu 
layers were relaxed. The PAW basis sets provide intrin- 
sically projections onto localized atomic orbitals, which 
we used to extract the hybridization functions (Eq. flS))) 
from our DFT calculations (for details see Refs. 12311241 ). 
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FIG. 2: (Color online) Hybridization functions and crystal 
fields for Co in bulk Cu (a) and as ad-atom on Cu (111) (b). 
In the upper panel the dynamical crystal field ed + Re A is 
shown. Lower panel: ImA. 



B. Impurity solver 

The impurity model ([3]) can be solved without approx- 
imations using continuous-time quantum Monte Carlo 
(CTQMC) algorithms. Both the interaction expansion 25 
and hybridization expansio n 26 ' 27 algorithms can treat 
multi-orbital systems with general four-fcrmion interac- 
tion terms. The hybridization expansion is advantageous 
in the case of a strongly interacting five orbital model, 
because the expansion in the hybridization leads to much 
lower perturbation orders than an expansion in the vari- 
ous interaction terms, and thus to much reduced compu- 
tational demands. For the interaction parameters of the 
Co impurities studied in this paper the order in the hy- 
bridization expansion was found to be a factor 20 lower 
compared to the interaction expansion. 17 Furthermore, if 
the hybridization function is diagonal in the orbital in- 
dices, as it is to a good approximation the case in the 
Co/Cu systems studied here, no sign problem appears. 
This is in contrast to the hybridization expansion, where 
correlated hopping terms were found to lead to severe 
sign cancellations^ 

Despite these advantages, a strong-coupling CTQMC 
simulation of a general five orbital model is still computa- 
tionally expensive, because the imaginary time evolution 
of the local impurity Hamiltonian H\ oc has to be com- 
puted exactly. The weight of a Monte Carlo configura- 
tion corresponding to In hybridization events (impurity 
creation and annihilation operators d\ and d a ) at times 
Ti < . . . < tL and 7"i < • • • < r n contains a factor— 
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where the trace and the operator products operate on a 
Hilbert space of H\ oc of dimension 1024. One strategy 
to evaluate this trace is to use the eigenbasis of H\ oc , 
in which the time-evolution operators are diagonal, and 
to order the eigenstates according to conserved quan- 
tum numbers. The operators da then acquire a block 
structure 2 ^ which reduces the effort for matrix-matrix 
multiplications. 



An alternative approach, which was found to be more 
efficient in the case of five orbital systems and low tem- 
perature, is the Krylov implementation. 18 In this algo- 
rithm one evaluates the trace factor in the occupation 
number basis. In this basis, the operators da are sparse 
and can easily be applied to any state, while the time- 
evolution operators become non-trivial dense matrices. 
However, we never evaluate the exponential of -ffi oc , but 
only exp(— H\ oc t)\v), for a given state \v). This can 
be done with only a small number of sparse matrix- 
vector multiplications, by using efficient Krylov space 
techniques. Details of the algorithm can be found in 
Ref. [lg. The computational advantage of this Krylov- 
space approach becomes particularly evident at low tem- 
peratures, where it is sufficient to consider the contri- 
butions of the low energy states in Tri oc [. . .], since the 
probability of the system to be in one of these states is 
very large (during the time evolution, all excited states 
are accessible) . If the outer trace is truncated in this way, 
the observables must be measured in the middle of the 
imaginary-time interval, i.e. at r = f3/2. 
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Co on Cu 
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7.3 
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0.96 
0.51 


7.6 
7.6 


1.00 
0.53 


7.4 
7.3 


0.96 
0.90 


7.7 
7.5 


1.00 
0.93 



TABLE I: Occupancies n and impurity spins S as obtained 
from our GGA and GGA+U calculations. Values obtained 
directly from the PAW projectors (n, S) and normalized by 
the integrated total Co d-electron DOS, J\f = J v(E)AE, are 
shown (n = n/Af, S z = S z /N). 



IV. RESULTS 

The DFT calculations yield the orbital dependent hy- 
bridization functions shown in Fig. [5] In bulk Cu, the 
environment of the Co impurities is cubically symmetric 
and the hybridization function decomposes into three- 
fold degenerate t 2g and twofold degenerate e g blocks. 
In the bulk symmetry forbids off-diagonal elements in 
the hybridization function. On the surface, the symme- 
try is reduced to C 3v . For Co on Cu the hybridization 
function decomposes into two twofold degenerate blocks 
transforming according to the E-irreducible representa- 
tion of C 3v (Ei (d xz , d yz ) and E 2 (d x 2_ y 2, d xy )) and the 
d z 2 -orbital transforming according to the A\ representa- 
tion. For Co on Cu the hybridization functions contain 
small off-diagonal matrix elements. These off-diagonal el- 
ements, which are smaller than the diagonal ones, will be 
neglected in our simulations. As a general trend, one can 
already see that the hybridization of the Co d-electrons 
is about twice larger in the bulk than on the surface. 

DFT calculations are also used to calculate the occu- 
pancy of the Co 3d impurity orbitals. To this end, we 
performed spin-polarized DFT calculations using GGA 
as well as GGA+U of Co in and on Cu with the full inter- 
action vertex defined via the average screened Coulomb 
interaction U = 4eV and the exchange parameter J = 
0.9 eV. We obtained the occupancies of the Co 3d orbitals 
derived from the PAW projectors n = n^ + nj_ and the 
impurity spin S z = h(rif — n^) and present them in Tab. 
HJ In all cases the average Co 3d occupancy suggested by 
our DFT calculations is between n — 1 and n — 8. For 
Co on Cu, the impurity spin is S z ~ 1 which is well in 
line with a d 8 configuration of the Co. In the bulk, the 
Co spin is S z « 1 in GGA+U and S z w 1/2 in GGA. 

In the following we study Co in and on Cu in the five- 
orbital Anderson impurity model formulation (Eq. ([I}). 
In this framework, the chemical potential has to be cho- 
sen to fix the occupancy of the Co d-orbitals. Due to the 
well know double-counting problem in LDA+DMFT type 
approaches,-- the precise chemical potential /i and the Co 
d-occupancy are not known. Therefore, we computed re- 
sults in a range of chemical potential values which yield 
a total d occupancy consistent with the estimates of the 
DFT calculations. For both systems the results of the 
DFT calculations predict a total density n < 8 and sug- 
gest a spin S~lor slightly below in the case of Co in Cu. 
For fi = 26, 27, 28 eV (Co in Cu) and n = 27, 28, 29 eV 
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FIG. 3: (Color online) Orbitally resolved DOS of the Co im- 
purities in bulk Cu (top) and on Co (111) (bottom) obtained 
from our QMC simulations at temperature T — 0.025 eV and 
chemical potential \i = 27 +V and [i = 28 eV, respectively. 



(Co on Cu) we obtain total densities and spins close to 
these DFT estimates. The values of both observables 
for the lowest simulation temperature T = 0.025 eV are 
presented in Table HU 



TABLE II: Total density and spin. 



System 


/i(eV) 


<n> 


(S) 


Co in Cu 


26 


7.51 + 0.07 


1.02 + 0.02 


Co in Cu 


27 


7.78 + 0.05 


0.92 ± 0.02 


Co in Cu 


28 


8.06 + 0.03 


0.817 ±0.007 


Co on Cu 


27 


7.76 + 0.05 


1.07 + 0.01 


Co on Cu 


28 


7.93 ± 0.05 


0.99 + 0.01 


Co on Cu 


29 


8.21 + 0.03 


0.860 + 0.007 
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FIG. 4: (Color online) DOS of the Co impurities in bulk Cu 
(top) and on Cu (111) (bottom) obtained from DFT (GGA 
and GGA+U) as well as QMC simulations at temperature 
T — 0.025 eV. QMC results obtained at different chemical 
potentials fi are shown. 



A. Quasi-particle spectra 

We now analyze the excitation spectra of the Co im- 
purities in order to understand the dominant physics at 
different energy scales. For a first, qualitative insight into 
the strength of many body renormalizations, we com- 
pare in Fig. |4] the Co 3d-electron DOS obtained from 
our DFT calculations to the Co 3d spectral functions ob- 
tained from analytical continuation of our QMC results. 

The non-spinpolarized GGA calculations used to de- 
termine the hybridization functions yield — by defini- 
tion — the LDOS corresponding to the Anderson model 
without two-particle interactions {U = J = OeV). For 
both, Co in and on Cu, the GGA DOS exhibits a peak 
near the Fermi level (E F = 0). The QMC DOS qualita- 
tively reproduces the GGA DOS for the case of Co in Cu. 



Here, the main difference between the two approaches is 
that QMC yields a peak near the Fermi level which is 
approximately twice narrower and shifted towards Ep. 
GGA+U accounts for the local Coulomb interactions at 
the Co atoms on a Hartree Fock level which leads to the 
destruction of the quasi-particle peak near Ep with all 
the spectral weight shifted to broad Hubbard bands. The 
comparison to the QMC results shows that this destruc- 
tion of the quasi-particle peak is unphysical. 

For Co on Cu the hybridization is weaker and the DOS 
from the QMC simulations exhibits both quasi-particle 
peaks near Ep as well as Hubbard type bands at higher 
energies. The reduction of spectral weight of the quasi- 
particle peak as compared to GGA is stronger here. 

The orbitally resolved DOS of Co in and on Co is shown 
in Fig. [3] For Co in Cu the DOS of the e g and the t^g 
orbitals is very similar particularly regarding the quasi- 
particle peak — despite the (energy dependent) crystal 
field splitting on the order of some 0.1 eV. 

The DOS of Co on Cu exhibits stronger differences 
between the Ei, E2, and A\ orbitals. The E2 orbitals, 
which spread out perpendicular to the z-axis, show the 
weakest hybridization effects, but even here, a quasi- 
particle peak appears in all orbitals. The appearance 
of low energy quasi-particle peaks in all orbitals is differ- 
ent from the behavior expected for a spin-1 two-channel 
Kondo model, where a low energy quasi-particle res- 
onance would be observed in two orbitals (four spin- 
orbitals) only. 

The DOS as obtained from our QMC calculations sug- 
gests a low-temperature Fermi liquid state involving all 
orbitals for both, Co in and on Cu. We investigate the 
nature of this state in the following sections by analyzing 
the self-energies obtained from QMC and the statistics 
of relevant atomic states. 



B. Low energy Fermi liquid 

If a Fermi liquid develops, the self-energy takes the 
form 



E(T, u) = X(T, 0) + £'(T, 0)w + 0(w 2 ) 



(7) 



with S(T, 0) and the first energy derivative X'(T, 0) being 
real for T — > 0. In this regime, the spectral weight Z 
associated with the quasi-particle peak is determined by 



Z = (1 - ReS'(O) - ReA'(0))" 



(8) 



Our QMC calculations yield the self-energy on the Mat- 
subara axis. Analytic continuation ui — > iui n shows that 
Fermi liquid behavior manifests itself on the Matsubara 
axis by 

ImS(T, iun) ~ ImS(T, 0) - Im£'(T, 0)w„ (9) 

at low frequencies with ImX(T, 0) ~ T 2 . 

We now compare these relations to the frequency and 
the temperature dependence of ImX(T, iu> n ) obtained 



from our QMC calculations. At the lowest accessible 
temperature, T = 0.025 eV, we obtained the Matsubara 
self-energies depicted for Co in and on Cu in Fig. [5] 

In both systems |ImS| clearly decreases as uj n — > 0, 
for all orbitals except for the Ei orbitals of Co on Cu 
at /J, = 27 eV. This is clearly different from the diverg- 
ing E(iw„) ~ -ri-, expected for the localized moment 
of an isolated atom. For Co in bulk Cu, the e g and 
tig orbitals exhibit very similar self-energies, whose low- 
energy behavior is consistent with the form expected for 
a Fermi-liquid (Eq. |§J). For Co on Cu, the self-energies 
differ considerably between the different orbitals with the 
E\ orbitals being least correlated and the Ei orbitals ex- 
hibiting the largest self-energies at low frequencies. Our 
results indicate that a Fermi liquid develops in all Co 
orbitals, also here, although the Kondo temperature ap- 
pears to be orbital dependent. 



C. Estimation of T K from QMC 

To define a Kondo temperature scale even in cases 
without well defined local moment at intermediate tem- 
peratures, we define Tr- through the width of the quasi- 
particle resonance in the single particle spectral function 
near Ep, which is measured in STM experiments. 

In our QMC simulations we determine Tk from the 
quasi-particle weight Z. The simulations yield the self- 
energy at the Matsubara frequencies u> n — n ~t . An- 
alytical continuation of Eq. © yields 



dImT,(iuj n ) 



diuj 71 



ReA'(O) 



(10) 
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and we use Eq. © to evaluate the derivative. In Fig. [5] 
we show the quasi-particle weight of Co in Cu and Co 
on Cu for the different types of orbitals as a function 
of the chemical potential jjl. The values of degenerate 
orbitals agree within an accuracy of 10~ 2 to which pre- 
cision we are listing them in Tab. IIIII for T = 0.025 eV. 
The systems whose spin is closest to S = 1 are found to 
have the lowest values of Z. Co in Cu clearly has higher 
quasi-particle weights compared to Co on Cu. As we 
will see, this results in a higher Kondo temperature Tk, 
which is also confirmed experimentally. In experiments 
using STM measurements the Kondo temperature has 
been found to be T K = 655^ ± 155AT = 0.056 ± 0.013 eV 
for Co in Cu 9 . and T K « 54 ± 5K = 0.0046 ± 0.0005 eV 
for Co on Cui&2&±5 

Following Hewson's derivation of a renormalized per- 
turbation theory of the Anderson models we use as a 
definition for the Kondo temperature: 



T K = ZlmA(O). 

4 



(11) 



The values computed according to Eq. (ITTj) are listed 
in Tab. IIIII for all impurity energy levels E a + /j, at tem- 
perature T — 0.025 eV. As for the quasi-particle weights 



we find the lowest values of Tjc for /i = 26 eV (Co in 
Cu) and fi = 28 eV (Co on Cu). Averaging over or- 
bitals we obtain T K = 0.118 eV (Co in Cu, fji = 26 eV, 
T = 0.025 eV) and T K = 0.016 eV (Co on Cu, n = 28 eV, 
T = 0.025 eV) with a ratio of T^ /T§ N = 7.4. This 
large difference between the two Kondo temperatures 
is in fair agreement with experiments, where a ratio 
of T™/T° N = 12.1 has been found. 9 As discussed in 
Section [V] the physical quantities which determine the 
Kondo temperature scale enter in the argument of an ex- 
ponential function, suggesting that a comparison of the 
logarithms log(T oxp )/log(TR-) is more appropriate when 
judging the predictive power of our first-principles sim- 
ulations. We find log(T cxp )/log(T K ) = 1.4 (Co in Cu, 
fi = 26 eV, T = 0.025 eV) and log(T cxp )/log(T K ) = 1.3 
(Co on Cu, n = 28 eV, T = 0.025 eV). 

For both systems our computed Kondo temperatures 
are higher than the experimentally determined values. 
This can be either due to the neglect of spin-orbit cou- 
pling effects which lifts degeneracies and narrows the low 
energy resonances or due to the Coulomb interactions 
being larger than the U = 4eV assumed here. 

The temperature dependence of the self-energy allows 
for an alternative test whether and at which energy scale 
a Fermi liquid emerges. This is illustrated for Co in and 
on Cu in Fig. [7] We find an almost temperature indepen- 
dent behavior of Im £(T, iw n ) for Co in Cu if T < 0.05 eV, 
which provides an estimate of Tk ~ 0.05 eV. In the case 
of Co on Cu Im£(T, iui n ) still evolves as one lowers the 
temperature from T = 0.05 eV to T = 0.025 eV, which 
suggests a lower Kondo temperature. 

According to Eq. (J9j) , the linear extrapolation of 
Im£(T, iuS) to iu! n — > yields the inverse lifetime Ht^ 1 = 
Im£(T, 0) of quasi-particles at the Fermi level if a Fermi 
liquid is formed. In Fig. [5] we plot this extrapolation for 



TABLE III: Kondo temperatures Tk computed from the 
quasi-particle weight Z and the hybridization function A(w) 
according to Eq. (flip at the lowest simulation temperature 
T = 0.025 eV. The experimental Kondo temperatures are 
T K = 0.056±0.013 eV (Co in Cu) and T K = 0.0046±0.0002 eV 
(Co on Cu)M 



System 


M(eV) 


E t 


-Im(A(0)) (eV) 


Z 


T K (cV) 


Co in Cu 
Co in Cu 


26 
26 


*2 S 

e 9 


0.43 ±0.01 
0.340 ± 0.009 


0.38 
0.39 


0.13 
0.10 


Co in Cu 
Co in Cu 


27 
27 


*2 S 

e 9 


0.43 ±0.01 
0.340 ± 0.009 


0.42 
0.47 


0.14 
0.12 


Co in Cu 
Co in Cu 


28 

28 


*2 9 
e 9 


0.43 ±0.01 
0.340 ± 0.009 


0.48 
0.56 


0.16 
0.15 



Co on Cu 


27 


E 2 


0.124 ±0.002 


0.06 


0.006 


Co on Cu 


27 


Ei 


0.226 ± 0.002 


0.19 


0.03 


Co on Cu 


27 


A! 


0.197 ±0.001 


0.08 


0.01 


Co on Cu 


28 


E 2 


0.124 ±0.002 


0.05 


0.005 


Co on Cu 


28 


Ei 


0.226 ±0.002 


0.15 


0.03 


Co on Cu 


28 


Ai 


0.197 ±0.001 


0.10 


0.01 


Co on Cu 


29 


E 2 


0.124 ±0.002 


0.14 


0.01 


Co on Cu 


29 


Ei 


0.226 ± 0.002 


0.26 


0.05 


Co on Cu 


29 


Ai 


0.197 ±0.001 


0.27 


0.04 
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FIG. 5: Orbitally resolved self-energies for Co in Cu (upper panel) and Co on Cu (lower panel). From left to right , ImE(iw) 
is shown for the chemical potentials fj, = 26, 27, 28 eV (Co in Cu) and /i = 27, 28, 29 eV (Co on Cu) . 
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FIG. 6: (Color online) Orbitally resolved quasi-particle weight 
of Co in Cu and Co on Cu at temperature T = 0.025 eV. 



Co in Cu, [i = 26 eV and show ImI](T, 0) as a function 
of temperature. We find a Im£(T, 0) ~ T 2 -behavior for 
both sets of orbitals, which corroborates the formation 
of a Fermi-liquid in the e g and t<i g orbitals at the lowest 
accessible temperatures. 



V. DISCUSSION 

Our QMC calculations showed that, for both Co in 
and on Cu Fermi liquids involving all orbitals of the Co 
impurity form at low temperatures. We now want to un- 
derstand these results on the basis of scaling arguments 
starting with (higher energy) charge fluctuations going 
to (lower energy) spin and orbital fluctuations. 



A. Charge fluctuations 

With our values of the Coulomb interaction strength 
in the local Hamiltonian (U = 4eV; J = 0.9 eV) the en- 
ergies of removing (E-) or adding (E + ) an electron to 
the impurity are £L ~ E + w 2 eV. A first order expan- 
sion in the hybridization gives a qualitative estimate of 
the role of charge fluctuations: 10 The norm of the admix- 
tures of dJ and d 9 configurations to a predominantly d 8, 
ground state of the impurity is approximately Af n ^& = 
-ilmA(0)j^. For Co in Cu, -ImA(O) w 0.4eV leads 
to N n =jt$, ~ 0.6. The hybridization of Co on Cu is about 
twice smaller, — Im A(0) « 0.2 eV, yielding a correspond- 
ingly smaller weight of non-d 8 configurations M n ^% ~ 0.3. 

Our QMC calculations allow to quantitatively mea- 
sure the charge fluctuations. In Fig. [9] we plot the "va- 
lence histogram"— for Co in/on Cu for different choices 
of the chemical potential and compare it to the "va- 
lence histogram" of the Slater determinant built from 
the lowest GGA eigenstates. The histogram shows the 
weights which the eigenstates in the different charge sec- 
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FIG. 7: (Color online) Upper panel: Im(E) of Co in Cu for /i = 28 eV at orbital energies Ei = — 0.288 eV {d xy ,d y z,d xz ) 
and Ei — — 0.44eV (d z 2,d x 2_ y 2) (from left to right). Lower panel: Im(E) of Co on Cu for /i = 29 eV at orbital energies 
Ei = -0.12 eV {d xy ,d x 2_ y 2), Ei = -0.39 eV (d yz ,d xz ) and Ei = -0.34 eV (d x i) (from left to right). 
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FIG. 8: (Color online) Temperature dependent quasi-particle 
lifetimes for Co in Cu at fi — 26 eV for the t^g (left) and e g 
orbitals (right). 



(k, 0.5 for Co in Cu and ss 0.3 for Co on Cu) agree sur- 
prisingly well with the simple perturbative estimate. 

A Schrieffer- Wolff decoupling of ionized impurity 
states to discuss the low energy physics can be justified if 
the weight of non-ci 8 configurations is N n =£% <C 1- Hence, 
an intermediate energy scale of well formed (unscreened) 
fluctuating spin or orbital moments at frozen impurity 
valency might be defined in the case of Co on Cu (111) 
but clearly not in the case of Co in Cu. 



B. Spin and orbital fluctuations 



tors n = 0, 1, ... 10 contribute to the partition function 
(via the outer trace in Eq. ([6])). In all cases the local 
Coulomb interaction in the QMC simulation leads to a 
narrower distribution of the occupancies as compared to 
the GGA valence histograms. This effect is most pro- 
nounced in the case of Co on Cu (111), where the d a 
configuration clearly dominates over the d 7 and d 9 con- 
figurations. For Co in Cu, there are still noticeable cor- 
relations and the narrowing of the valance histogram as 
compared to the GGA case but the d s configuration con- 
tributes only about 50% in the QMC simulations, with 
significant weight coming from the d 7 , d 9 and even the d 6 
and d 10 configurations. The measured values for Af n =£8 



While for Co in Cu charge, spin and orbital fluctu- 
ations will be present down to lowest energies, for Co 
on Cu only fluctuations of the orbital and the spin de- 
gree of freedom are expected to dominate the low energy 
physics. To further investigate to which extent orbital 
and spin fluctuations might determine the low energy be- 
havior of the impurity we estimate Kondo temperatures 
within simplified models and compare these estimates to 
our QMC calculations as well as experiments. 

Assuming well-defined ma gnet ic moments (Af n ^s "C 1) 
a scaling analysis (cf. Ref. llOh allows to estimate the 
Kondo scale analytically in simplified situations: If nei- 
ther Hund's rule coupling nor crystal field splitting or any 
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other symmetry breaking terms were present the spin and 
the orbital degrees of freedom of the Co impurities could 
fluctuate freely and independently. This would lead to a 
Kondo temperature^ Tk ~ A) exp [— \/2N n ^&\, where 
D a ~ min(E±, A) is related to the impurity charging en- 
ergies (E±) and the electronic bandwidth A. Dq is on 
the order of several eV. With Dq = U/2 — 2 eV we es- 
timate Tk w 0.4_D = 0.9 eV in the case of Co in Cu 
and T K « 0.2£> = 0.4 eV in the case of Co on Cu. This 
is in both cases (at least) an order of magnitude larger 
than the Kondo temperatures obtained from QMC and 
the experimentally measured Kondo temperatures. 

The opposite limit is given by the case with strong 
Hund's rule coupling and supressed orbital fluctuations. 
Without orbital fluctuations, but still disregarding the 
Hund's coupling J, the Kondo temperature readsi^ 



T 



K 



Dq exp 



ttU 



= Dq exp 



5 



which 



8lmA(o) 

would lead to T K = 0.02D « 0.04 eV for Co in Cu. The 
Hund's rule coupling reduces the Kondo temperature 12 
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FIG. 10: (Color online) Comparison of the DOS (left) and 
quasi-particle weight Z (right) of Co on Cu at U — 4eV for 
the two values J = 0.0 eV and J = 0.9 eV. 



further to T K = T K (T K / JhS) 23 - 1 . With S = 1 and 
Jh = 0.9 eV this would lead to T K » 0.002 eV. For 
Co on Cu, the assumption of an orbital singlet yields 
T K = 0.0004D s* 0.001 eV and Hund's rule coupling fur- 
ther reduces the Kondo temperature to T K « 1 fieV. This 
limit thus yields Kondo temperatures which are orders of 
magnitude smaller than those obtained in our QMC cal- 
culations as well as the Kondo temperatures measured 
experimentally for Co in and on Cu. 

It is thus the successive locking of the impurity elec- 
trons to a large spin by the Hund's rule coupling and 
the partial freezing out of orbital fluctuations that deter- 
mines the on-set of Fermi liquid behavior and the Kondo 
temperature in realistic systems like Co in or on Cu. 

With this in mind it is instructive to analyze the in- 
fluence of a static crystal field on the energy spectrum 
of otherwise isolated Co atoms. Without crystal fields, 
in a d 8 configuration our local Coulomb interaction (U — 
4eV; J = 0.9 eV) yields a 21 fold degenerate L = 3, S = 1 
ground state which is separated from the L = 2, S = 
multiplet by an energy of El=2,s=o — 1-3 eV. This is 
clearly larger than the crystal field acting on the Co im- 
purities (Fig. [2]): The cubic crystal field (evaluated at 
the Fermi level) of Co in Cu leads to the e g states be- 
ing 0.18 eV higher in energy than the ti g states. In this 
crystal field, the resulting d 8 ground state is an orbital 
singlet. Excitations to higher crystal field split states re- 
quire energies on the order of 0.2 eV. This is larger but 
comparable in order of magnitude to the Kondo temper- 
atures obtained in experiments and simulation. However, 
fluctuations to these higher crystal field split states must 
be taken into account to explain the characteristic tem- 
perature of the low energy Fermi liquid formed at Co 
impurities in Cu. 

In our model of Co on Cu, the static crystal fields 
also lift the degeneracy of the ground state multiplet 
but a double degeneracy in the orbital space remains. 
Excitations to higher crystal field split states require 
0.03 — 0.08 eV. In this model, even the ground state mul- 
tiplet allows for fluctuations of the orbital degree of free- 
dom. 

In order to examine the effect of constraining orbital 
fluctuations we consider the case of Co on Cu which 
exhibits a strong reduction of the quasi-particle peak 
compared to the GGA spectral function representing the 
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U = 0, J = case. Turning off the Hund's coupling J al- 
lows the orbital and spin degrees of freedom to fluctuate 
more freely and given the scaling considerations should 
result in a higher Kondo temperature as well as a broader 
quasi-particle peak. We study the effect of J = for Co 
on Cu, T — 0.025 eV, /j, = 29 eV and present the com- 
parison of the quasi-particle weight and peak in Fig. 1101 
In line with our statement that the Kondo temperature 
is determined by the locking of the impurity electrons to 
a larger spin and possible restrictions of the orbital fluc- 
tuations, we find a broadening of the quasi-particle peak 
and increase of the quasi-particle weight Z as J — > 0. 



understanding of magnetic nanostructures based on 3d 
adatoms on surfaces as well as magnetic impurities in 
bulk metals require us to account for the orbital degrees 
of freedom. Dynamical mean field theory provides a link 
between quantum impurity problems and extended lat- 
tices of atoms which are subject to strong electron cor- 
relations. 

It remains thus a future challenge to understand how 
the orbital degree of freedom controls the quenching of 
magnetic moments and eventually the formation of low 
energy Fermi liquids in realistic extended correlated elec- 
tron systems as well as magnetic nanostructures. 



VI. CONCLUSIONS 

For Co in and on Cu we found that a Fermi liquid is 
formed at low T involving all impurity d-orbitals. The 
example of Co on Cu shows that the characteristic tem- 
perature, Tk, associated with the onset of Fermi liquid 
behavior can differ between the impurity orbitals. The 
comparison of our QMC calculations and scaling argu- 
ments further demonstrates that fluctuations in the or- 
bital degree of freedom and Hund's rule coupling are cru- 
cial in determining Tk in realistic systems. This is well 
beyond the physics of simple "spin-only" models. The 
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